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Whilst silicon is a promising material for quantum computation, the degeneracy of the conduction 
band minima (valleys) must be lifted with a splitting sufficient to ensure formation of well-defined 
and long-lived spin qubits. Here we demonstrate that valley separation can be accurately tuned via 
electrostatic gate control in a metal-oxide-semiconductor quantum dot, providing splittings spanning 
0.3-0.8 meV. The splitting varies linearly with applied electric field, with a ratio in agreement with 
atomistic tight-binding predictions. We demonstrate single-shot spin readout and measure the 
spin relaxation for different valley configurations and dot occupancies, finding one-electron lifetimes 
exceeding 2 seconds. Spin relaxation occurs via phonon emission due to spin-orbit coupling between 
the valley states, a process not previously anticipated for silicon quantum dots. An analytical theory 
describes the magnetic field dependence of the relaxation rate, including the presence of a dramatic 
rate enhancement (or hot-spot) when Zeeman and valley splittings coincide. 



I. INTRODUCTION 

Silicon is at the heart of all modern microelectronics. 
Its properties have allowed the semiconductor industry 
to follow Moore's law for nearly half a century, delivering 
nowadays billions of nanometre-scale transistors per chip. 
Remarkably, silicon is also an ideal material to manipu- 
late quantum information encoded in individual electron 
spin4^^. This is a consequence of the weak spin-orbit 
coupling and the existence of an abundant spin-zero iso- 
tope, which can be further enriched to obtain a "semicon- 
ductor vacuum" in which an electron spin can preserve a 
coherent quantum superposition state for exceptionally 
long times^. 

In order to define a robust spin-^/2 qubit Hilbert space, 
it is necessary that the energy scale of the two-level sys- 
tem is well separated from higher excitations. In this 
respect, a major challenge for the use of silicon is repre- 
sented by the multi- valley nature of its conduction band. 
In a bulk silicon crystal the conduction band minima are 
six- fold degenerate, but in a two-dimensional electron gas 
(2DEG), the degeneracy is broken^ into a two- fold degen- 
erate ground state (F valleys) and a four-fold degenerate 
excited state (A valleys), due to vertical confinement of 
electrons with different effective mass along the longitudi- 
nal and transverse directions, respectively. Furthermore, 
the F valley degeneracy is generally lifted by a sharp 
perpendicular potentiaP^ and the relevant energy sepa- 
ration is termed the valley splitting (VS). 

The valley splitting depends on physics at the atomic 
scald^^HI^ (e.g. roughness, alloy and interface disorder), 
and so it is not surprising that experiments have revealed 
a large variability of spli ttings among devices, ranging 
from hundreds of jieV^^^^^ up to tens of meV in excep- 
tional cases^^. At present, the lack of a reliable experi- 



mental strategy to achieve control over the VS is driving 
an intense research effort for the development of devices 
that can assure robust electron spin q ubits by minimising 
multi-valley detrimental effecta^^^, or even exploit the 
valley degree of freedomP^^ for new types of qubits. 

Another crucial parameter to assess the suitability of a 
physical system to encode spin-based qubits is the relax- 
ation time of spin excited states (Ti). Spin lifetimes have 
been measur ed for gate-defined Si quantum dots (QDs)P^, 
Si/SiGe Qp j^^ l ^^ j donors in Si^^, reporting values 
which span from a few milliseconds to a few seconds. 
Furthermore, the dependence of the spin relaxation rate 
(Tj~^) on an externally- applied magnetic field {B) has 
been investigated. Different mechanisms apply to donors 
and QDs, accounting for observed T^^ ex and B^ de- 
pendencies^, respectively. In principle, T^^{B) depends 
on the valley configuration and the details of the excited 
states above the spin ground state. However, until now, 
no experimental observation of the effects of a variable 
VS on the relaxation rate has been reported. 

Here, we demonstrate for the first time that the valley 
splitting in a silicon QD can be finely tuned by direct 
control of an electrostatic gate potential. We find that 
the dependence of the VS on vertical electric field at the 
Si/Si02 interface is strikingly linear, and show that its 
tunability is in excellent agreement with atomistic tight- 
binding predictions. We demonstrate accurate control of 
the VS over a range of about 500 /ieV and use it to ex- 
plore the physics of spin relaxation for different QD occu- 
pancies {N=l, 2, 3). We probe both the regime where the 
VS is much larger than the Zeeman splitting at all mag- 
netic fields and that where the valley and spin splittings 
are comparable. We observe a dramatic enhancement of 
the spin decay rate (relaxation hot-spot) when spin and 
valley splittings coincide. To our knowledge, such hot- 
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FIG. 1: Device architecture and addition energy spectrum. 

(a) Schematic (top view) of the device's gate layout. Differ- 
ent colours represent different layers within the gate stack. 

(b) Schematic diagram of the single-lead QD (left) and SET 
detector (right). Regions where an electron layer is formed are 
coloured in orange. The readout signal (/set) is sensitive to 
the QD charge state due to the QD/SET capacitive coupling 
(C'cpi). (c) Device cross-sectional schematic. An electron layer 
is formed underneath the positively biased gates: Rl and R2 
define the QD reservoir; P controls the QD population; and 
ST the sensor's island. The Si02 layer (in purple) thickness 
and plunger gate width are indicated, (d) Energy diagram 
showing qualitatively the conduction band profile in the de- 
vice. Electrons accumulate wherever the gate bias lowers the 
conduction band below the Fermi level, Ef. (e) Charging en- 
ergy as a function of electron number. Spikes corresponding 
to complete 2D shell filling are observed, (f) Schematic of elec- 
tron filling for two- valley 2D Fock-Darwin states. Each state 
can hold two electrons of antiparallel spin and is identified 
by a pair of quantum numbers (n, /) and its valley occupancy 
(v). 



spots (HS) have been predicted for relaxation involving 
orbital states^^ but they are yet to be observed. We 
develop an analytic theory that explains the 5-field de- 
pendence of the relaxation rates and the details of the 
relaxation HS in terms of admixing of spin- valley states. 
This mechanism is seen to be significantly more promi- 
nent than the conventional spin-orbit hybridization^. 



II. QD ADDITION SPECTRUM 

Our device is fabricated using a multi-level gated 
metal-oxide-semiconductor (MOS) technology'^, and its 
architecture is depicted in Fig. [TJa)-(d). A quantum dot 
is formed under gate P by applying a positive bias voltage 
to induce an electron accumulation layer. Strong planar 



confinement for the dot's potential well is achieved by 
negatively biasing gates B, CI and C2. A 2DEG reser- 
voir is also induced by positively biasing gates Rl and R2, 
and the QD occupancy can be modified by inducing elec- 
trons to tunnel between this reservoir and the dot. The 
remaining gates, namely SBl, SB2, ST, are employed 
to define a single-electron transistor (SET), capacitively 
coupled to the QD and used as a read-out device. The 
high flexibility of our design would allow us to use the 
same device also as a (single-lead) double dot structure 
by rearranging the gate bias (e.g. dots can be formed 
under gates B and C2). However, in this work, we only 
present results relevant to the single-dot configuration. 

In order to characterize the addition spectrum of the 
QD, we make use of a technique previously developed 
for GaAs-based systems which combines charge detection 
and gate pulsing^^. There is no direct transport through 
the single- lead QD and, therefore, addition/removal of 
charge is only detected via modifications in the SET cur- 
rent. In particular, charge transitions are detected as 
current peaks in the SET signal whenever the QD en- 
ergy eigenstates come into resonance with the reservoir's 
Fermi level. Note that the SET-QD coupling is merely 
capacitive (via Ccpi) and electrons do not tunnel between 
them. In order to maximize charge sensitivity in the de- 
tector, we employ dynamic voltage compensation^^ on 
different gates which makes our read-out signal virtually 
unaffected by slow charge drifts and random charge re- 
arrangements. 

Fig. 1(e) illustrates the addition energy spectrum for 
the first fourteen electron additions to the QD. There 
is very little variation of charging energy {Eq) for high 
occupancies, i.e. Eq ~ 11 meV for > 9. However, 
by decreasing the electron number, the charging energy 
steadily increases, as expected when the dot size is signif- 
icantly affected by the electron number. This evidently 
indicates that the few-electron regime has been achieved. 
Most interestingly, the energy spectrum shows peaks for 
the addition of the 5th and 13th electrons. The extra 
addition energy needed for those transitions can be at- 
tributed to complete filling of the first and second or- 
bital shells. As illustrated in Fig. 1(f), this is consistent 
with the energy spectrum of two-valley 2D Fock-Darwin 
stated, where the first and second orbital shells hold 
4 and 8 electrons respectively. This confirms that we 
can probe the occupancy until the last electron. To our 
knowledge, such a neat manifestation of two-dimensional 
shell structure has been observed before only in InGaAs 
dots^^ and in Si/SiGe dots^^. It is remarkable that our 
MOS-based planar architecture can show such a degree 
of circular symmetry. This is likely due to both our gate 
design (leading to symmetric dot confinement) and ex- 
cellent quality of the thin gate oxide. 
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FIG. 2: Spin readout and relaxation rates for single-electron occupancy, (a) Schematic diagram showing the effect of the 
3-level pulse sequence on the electro-chemical potential of the dot. Energy levels in the QD are Zeeman split according to spin 
polarisation and valley degeneracy is lifted. For clarity, only lower valley states are shown to be loaded/unloaded, (b) Pulsing 
sequence (top) for the single-shot spin readout and normalised SET signal for spin-up (middle) and spin-down (bottom), (c) 
Energy diagram of the 1-electron spin- valley states as a function of B-field. Maximum mixing of spin and valley degrees 
of freedom occurs at the anticrossing point where Zeeman and valley splittings coincide. Relevant relaxation processes are 
sketched, (d) Relaxation rates as a function of magnetic field for different valley splittings. Data points for ^vs = 0.75meV, 
Eys = 0.33meV are shown as green and red circles, respectively. Dashed lines are the calculated relaxation rates fitted with 
r=1.7 nm (green), r=l.l nm (red). 



III. SPIN- VALLEY LIFETIMES 



In order to measure the spin state of individual elec- 
trons in the QD, we use an energy-selective readout tech- 
nique^^. The readout protocol consists of three phases 
clocked by a three-level pulsed voltage applied to gate 
P which directly controls the dot's electrochemical po- 
tential [see Fig. [2] (a)]. Firstly, an electron of unknown 
spin is loaded into the dot causing a sudden decrease 



in the sensor current. Next, the potential of the dot is 
lifted so that the Fermi level of the reservoir lies between 
the spin- up and spin-down states of the dot, meaning 
that a spin-up electron can tunnel off the dot while a 
spin-down electron is blocked. This is the read phase, 
during which the presence of a spin-up state would be 
signalled by a current transient (spin-up tunnels out and 
then spin-down tunnels in) whereas a spin-down electron 
would lead to no current modification. Finally, the dot's 
potential is further lifted to allow the electron to tunnel 
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off, regardless of its spin orientation. In Fig.[2](b) single- 
shot traces for both spin-up (in blue) and spin-down (in 
green) detection are plotted. The longer the system is 
held in the load phase before performing a read opera- 
tion, the more likely it is for the spin-up excited state to 
decay to the spin-down ground state. Thus, by varying 
the length of the load phase and monitoring the probabil- 
ity of detecting a spin-up electron, we can determine^^ 
the spin lifetime, Ti. In our experiments the 5-field is 
directed along the [110] crystallographic axis. As shown 
in Fig. [2](d), we observe a wide range of spin lifetimes as 
a function of magnetic field, with lifetimes as long as 2.6s 
at the lowest fields studied, 5=1. 25T. These are some of 
the largest lifetimes observed to date in silicon quantum 
dotPl. 

A key focus of our experiment was to electrostatically 
tune the valley energy separation and measure relaxation 
rates in different valley configurations and QD electron 
occupancies. As we show below, our data definitively in- 
dicate that excited valley states play a critical role in the 
spin relaxation processes. We develop a theory to explain 
how changes in the valley splitting affects the spin- valley 
state mixing and leads to the observed relaxation times. 

As we detail in Section \P7\ we have attained accurate 
gate control of the valley splitting, allowing us to tune 
it over a range of hundreds of /ieV. This permits us to 
conduct experiments in regimes where the valley split- 
ting {Eys) is either larger or smaller than the Zeeman 
spin splitting (£^z), depending on the magnitude of the 
magnetic field [see Fig. |2](c)]. 

Fig. [2] (d) presents measurements of spin relaxation 
rates as a function of magnetic field for two valley split- 
ting values at a fixed dot population of = 1. We start 
by examining a configuration where the valley separation 
is larger than the spin splitting at all fields (green data 
set). In other words, we operate in a regime for which 



Evs > Ez= giiBB 



(1) 



where g is the electron gyromagnetic ratio, /iB is the Bohr 
magneton and B is the applied in-plane magnetic field. 

For ^vs=0-75 meV [green data in Fig. [2] (d)], we ob- 
serve a monotonic increase of the rate with respect to B 
that becomes increasingly fast as Ez approaches £^vs • In 
our experimental conditions, the T^^ cx B^ depend ences 
for known mech anisms i n silicon should not appl}^^'^, 
while predicted^^SIEzEHl rates cx B^ do not explain the 
experimental data. 

By decreasing the valley separation to Esjs =0.33 meV, 
we can achieve the condition where the Zeeman splitting 
matches or exceeds the valley splitting. The red data 
in Fig. [2](d) illustrate the situation where inequality (1) 
only holds for B < 2.8T. When Ez = Eys (i.e. for 
5=2. 8T), a spike in the relaxation rate occurs. Relax- 
ation hot-spots have been previously predicted to occur 
for spin rela xation inv olving orbital states in single and 
coupled QDp^a^^. To our knowledge, this is the first 
experimental observation of such a phenomenon. 

In order to understand the relaxation mechanisms, we 



have developed a model that takes into account the per- 
turbations in pure spin states due to spin-orbit coupling 
(SOC), yielding eigenstates which are admixtures of spin 
and valley states. The four lowest spin- valley states [see 
Fig. [2|;c)] are defined as |1) = |2) = |v-,t), 

1 3) =]'U+,|), 1 4) = |'y+,t)- These states are considered 
to be only very weakly affected by higher excitations, 
such as orbital levels which are at least 8 meV above the 
ground state in our devic^^. In practice, mixing to a 2p- 
like orbital state leads to a B^ dependence in T^^ and is, 
therefore, important mainly for high 5-fields (above the 
anticrossing point). At lower fields, the prominent mech- 
anism is the spin-valley admixing, which we now discuss 
in detail. 

The relaxation between pure spin states is forbidden 
because the electron-phonon interaction does not involve 
spin fiipping. However, in the presence of interface disor- 
der, SOC can mix states that contain both the valley and 
spin degrees of freedom, thus permitting phonon-induced 
relaxation. Indeed, in the non-ideal case of QDs with a 
disordered interface, roughness can perturb the envelope 
function of both valleys (otherwise identical for ideal in- 
terfaces) and allows one to assume non-zero dipole matrix 

elements, such as r ^ = {v-\r\v-^)^ r = {v-\r\v-), 

r++ = {v-^\r\v-^). By means of perturbation theory, we 
define renormalized excited states |2) and |3) that relax 
to the ground state |1). The leading order corrections to 
the wavefunctions are given by: 



|2)(°^=sin||2)-cos||3) 



(2) 



|3) 



(0) 



cos- |2) +sin- |3) 
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(3) 



where cos^ 



[^]i/2, sin^ = [^]i/2^ and a 



—8/^8'^ + is an expression involving the detuning 
from the anticrossing point, 5 = Eys — Ez^ and the en- 
ergy splitting at the anticrossing, A^. 

By evaluating the relaxation rate via the electron- 
phonon deformation potential interaction, we obtain the 
rate below the anticrossing as: 



■ 21 



cos 



27 



(4) 



where Ty^y = r^/^(A£^^/^, r ^) is the pure valley relax- 
ation rate that depends on the energy difference relevant 
to the transition and the dipole matrix element. The ex- 
perimental condition for which the hot-spot occurs (i.e. 
Eys = Ez) is modelled as an anticrossing point for the 
mixed states |2) and |3). At that point, spin relaxation 
is maximized and approaches the valley relaxation 
rate, as (5 ^ in eq. Q. 

Above the anticrossing (i.e. Eys < Ez), the relevant 
relaxation transitions are |3) ^ |1) and |3) ^ |2) (the 
subsequent decay |2) |1) is in the form of a fast inter- 
valley transition and is, therefore, neglected). The ana- 
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FIG. 3: Spin- valley relaxation for 2-electron occupancy, (a) Energy diagrams of the 2-electron spin- valley states in the dot's 
potential well (left) and as a function of B-field (right). Dashed line indicates that T+ is not accessible (see text), (b) Relaxation 
rate as a function of 5-field for N — 2 and ^vs = 0.58meV. Red (blue) crosses represent data points at fields smaller (larger) 
than the anticrossing point. Dashed lines are the calculated rates fitted with r^^=4.76 nm. Inset: Data from the main graph 
re-plotted as a function of the modulus of the detuning energy, 5. Points of equal absolute detuning have nearly the same decay 
rates. 5^ dependence (grey line) is a guide for the eye. 



lytical formulations of these contributions read: 
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(5) 



2 'T 2 'T 

Tga = sin - cos - r„/„(Ai;„/„, r_ 



■r++) (6) 



The dashed lines in Fig. 2(d) show the calculated re- 
laxation rates relevant to the two experimental values of 
£^vs discussed. We use dipole matrix elements as a single 

free parameter by assuming \r ^| \r — r^+l = r. 

A least-square fit to the experimental data is performed 
by fixing the SOC strength to {jSn — oir)^ 45 — 60 m/s 
(justif ied b y the high electric field 3 x lO^V/m, see 
Refs. P^), where Pd {olr) is the Dresselhaus (Rashba) 
strength. The fit then extracts a dipole size r ^2-3 nm 
for all three values of £^vs- 



The good agreement between the calculations and ex- 
periment, as well as the presence of a hot-spot at the 
point of degeneracy between Zeeman and valley split- 
ting, provide strong evidence of our ansatz that the 
spin relaxation is predominantly due to a new mecha- 
nism: that of mixing with the excited valley states via 
Rashba/Dresselhaus-like SOC. 

We now examine the case N = 2 electrons, and inves- 
tigate the dependence of the relaxation rate on the mag- 
netic field at a fixed valley splitting (£^vs=0-58 meV). 
We note that the energy levels accessible for loading the 
second electron in the dot, when the N = 1 spin-down 
ground state is already occupied, are either the singlet 
{S) or the two lower triplets (T_,To), while the higher 
triplet (T+) would require a spin- flip and is, therefore, not 
readily accessible [see Fig. |3ja)]. In general, for triplet 
states, the antisymmetry of the two-electron wavefunc- 
tion requires one electron to occupy a higher energy state. 
For our multi- valley QD [Fig. [l] (f)], this requirement is 
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fulfilled when the two electrons occupy different valley 
states [see Fig. [3](a)]. For low fields, the ground state is 
S and the triplets have higher energies. This results in 
excited states (triplets) that extend over two valleys and 
relax to single- valley ground states (singlets). 

As the magnetic field is increased, S and T_ undergo 
an avoided crossing {B = Bst)^ and then T_ becomes the 
ground state. We adjust the levels of our pulsed readout 
protocol so that, during the load phase only S and T_ are 
below the reservoir's Fermi energy, while, during the read 
phase, the Fermi energy is positioned within the singlet- 
triplet (ST) energy gap. As a consequence, for B < Bst 
{B > Bst) a T_ (S) state would be signalled with a 
current transient, and relaxation rates can be extracted 
as for the N = 1 occupancy. The experimental relax- 
ation rates in Fig. [sj^b) show a strongly non-monotonic 
behaviour, approaching an absolute minimum at the an- 
ticrossing point (5st=5T). The trend is strikingly sym- 
metric, as can be appreciated when T^^ is plotted against 
the detuning energy, as shown in the inset. This symme- 
try is reflected in the QD energy spectrum [Fig. [sj^a)], as 
far as the S and T_ states are concerned. For B < Bst^ 
the ST energy gap decreases with increasing 5, resulting 
in slower relaxation rates. By contrast, for B > ^st, the 
ST energy gap increases for increasing field, and so does 
the relaxation rate. 

To model the 2-electron case, we build the wavefunc- 
tions for S and T_ from the single-particle states by 
considering the Coulomb interaction as a perturbing av- 
eraged field. The corresponding states are defined as 
|l)2e = \v-,V-,S), |2)2e = |v_,v+,T_). Next, the addi- 
tional perturbation given by SOC leads to renormalized 
eigenstates which are admixtures of singlet and triplet, 
and whose leading order approximations read: 

|T)^°j=sin||l)2,-cos||2),, (7) 

|2)^?=COs||l)2e + Sm||2)2e («) 

being similar forms to eq. ([2| ([3|. By evaluating 
the electron-phonon Hamiltonian matrix element for the 
transition between these states, one finds that it coincides 
with its 1-electron counterpart for |3)^^ |^)ie* There- 
fore, we can conclude that the corresponding relaxation 
rates have the same functional form as those derived in 
eq. ([6|, i.e. F^j = F^^, although the matrix elements for 
the two cases will be different. Dashed lines in Fig. [3](b) 
represent the calculated rates which are fitted to the ex- 
perimental data similarly to the case N = 1. Once again, 
the model convincingly reproduces the main features of 
the experimental trend, in particular the rates for fields 
away from the anticrossing together with the symmetry 
of the characteristics with respect to Bst- Further work 
may be needed to improve the fit in the vicinity of the 
anticrossing point. 

We also measured the relaxation rates for N = 3 
electrons (not shown) and, interestingly, observed nearly 



identical rates to the 1-electron case. This may be ex- 
plained in terms of a negligible effect of the Coulomb 
interaction on the valley spectrunPH^. 



IV. VALLEY SPLITTING CONTROL 

We now turn to the experimental demonstration of ac- 
curate control of the valley splitting, ^^vs, via electro- 
static gating. To determine £^vs, we use two different 
experimental approaches. One utilizes the rapid increase 
in spin relaxation at the hot-spot, and is applicable in 
the low magnetic field regime. The other is based on 
magnetospectroscopy, and is relevant for high fields. 

The first technique stems from the fact that the hot- 
spot (HS) can be reliably detected by monitoring the 
spin-up probability as a function of magnetic field. In 
Fig. [4] (a), we show measurements of the spin- up proba- 
bility performed with the same method as the one used 
to evaluate spin lifetimes. We see that the probability 
of detecting a spin-up electron decreseas significantly at 
some magnetic fields. A sudden drop of the spin- up frac- 
tion in a narrow range of field identifies the increase in 
relaxation rate associated with the HS. Given that val- 
ley and Zeeman splittings coincide at the HS, one can 
extract the valley separation as Eys = ^/^b^hs^ where 
Bus is defined as the field at which the HS is observed. 
For varying gate voltage configurations, we scan B in 
the range 2.8T< B <5T, and identify Byls by setting 
an arbitrary probability threshold [green shaded area in 
Fig. |4ja)] below which the HS is assumed to occur. The 
use of this technique is limited to B < 5T because the life- 
time drop at the HS can be therein confidently assessed. 
At higher fields the relaxation becomes increasingly fast 
and its enhancement at the HS is indistinguishable within 
our measurement bandwidth ( ^ lOKHz). 

In order to evaluate Eys at higher magnetic fields, 
we use a more conventional magneto-spectroscopic ap- 
proach, as shown in Fig. |4] (b). By employing the same 
gate-pulsed technique used for the charge stability ex- 
periments, we focus on the singlet-triplet ground state 
transition as we load the second electron into the dot 
(i.e. = 1 ^ 2 transition) in the range 5T< B <6.5T. 
This is clearly identified as the point where the S (light 
grey feature) and T_ (dark grey feature) states cross. 
Here, Eys = QIJ^bBst^ as seen in Fig. [sj^a). 

The data points in Fig. |4] (c) represent the measured 
valley separation as a function of T/p, obtained by means 
of the aforementioned techniques. The solid line fit shows 
remarkable consistency between the two sets of data and 
reveals that E-ys depends linearly on the gate voltage over 
a range of nearly 500 /ieV, with a slope of 640 /ieV/V. 
To our knowledge, this is the first demonstration of the 
ability to accurately tune the valley splitting electrostat- 
ically in a silicon device. 

A linear dependence of the valley splitting with re- 
spect to the vertical electric field has b een pr edicted for 
2DEG systems via effective mass theor30^tioi ^ similar 




FIG. 4: Tunability of the valley splitting via gate- voltage control, (a) Spin- up probability as a function of magnetic field for 
Vp = 1.59V. The occurrence of a HS minimizes the spin-up fraction and allows one to extract Bus (shaded area) and, in turn, 
EyS' (b) Pulsed- voltage magnetospectroscopy showing dls^T/dVp at A/" = 1 ^ 2 transition. A square pulse of amplitude 16 
mV at 287 Hz is applied to gate P. The evolution of the energy difference between the singlet state (light grey) and the triplet 
state (dark grey) allows one to extract Bst (dashed line) and, in turn, ^vs- (c) Valley splitting as a function of plunger gate 
voltage (bottom axis) and modulus of interface vertical electric field in the QD (top axis). Blue and red dots show the valley 
separation (left axis) measured with HS and magnetospectroscopic techniques, respectively. The linear fit (solid line) indicates 
a valley tunability of 0.640 meV/V. Dashed lines are calculations performed via atomistic (purple) and effective mass (green) 
methods (right axis). Bottom inset: Table comparing calculated and experimental tunabilities. Top inset: TCAD simulation 
of the conduction band profile at Vp = 1.4V in the (x, z) plane. Negative energy region in red reveals where the dot is formed. 



dependence for MOS-based QD&^^ has also been reported 
by employing atomistic tight-binding calculations^^. In 
order to compare our experimental finding with the the- 
oretical predictions, we simulate the vertical electric field 
(Fz) in the vicinity of the dot for the range of gate volt- 
ages used in the experiments. We employ the commercial 
semiconductor software ISE-TCAD^'' to model the device 
electrostatic potential, and thereby the electric fields in 
the nanostructure. For this purpose, TCAD solves the 
Poisson equation with an approximation of Newton's it- 
erative method^^ to obtain convergence at low tempera- 
tures. 

The spatial extent of the dot is identified by regions 
where the calculated conduction band energy drops be- 
low the Fermi level [red area in the top inset of Fig. 



4(c)]. Note that our calculations are performed on a 
three-dimensional geometry identical to the real device 
with the only free parameter being the amount of offset 
interface charge. This is adjusted to match the experi- 
mental threshold voltage of the device ( V|;h =0-625 V). 

The computed variation of interface electric field with 
gate voltage Vp is used to determine the valley splitting 
according to both the atomistic^^ and effective mas^ pre- 
dictions. Dashed lines in Fig. [4] (c) depict the trends for 
both approaches, with both exceeding by more than 1 
meV the measured values. Despite this offset, the atom- 
istic calculations give a tunability of the valley splitting 
with gate voltage, A£^vs/^^P7 in good agreement with 
the experiments. The calculated value of 597 /ieV/V 
agrees with the measured value to within less than 7%. 



8 



The value of 541 fieV/V calculated using the effective 
mass approach reveals a larger deviation 15%) from 
the experiments. The presence of an offset in the com- 
puted valley splitting may be due to the contribution of 
surface roughness that is not accounted for in the models, 
and is th ought to be responsible for a global reduction of 
^yj811Qi ri2.49,5Q. emphasize, however, that the gate 
tunability would remain robust against this effect, which 
is not dependent on electric field. 

V. CONCLUSIONS 

In this work, we have shown that the valley splitting in 
a silicon device can be electrostatically controlled by sim- 
ple tuning of the gate bias. We used this valley splitting 
control, together with spin relaxation measurements, to 
explore the interplay between spin and valley levels in a 
few-electron quantum dot. 

The relaxation rates for a one-electron system exhibit 
a dramatic hot-spot enhancement when the spin Zeeman 
energy equals the valley splitting, while for a two-electron 
system the rates reach a minimum at this condition. We 
found that the known mechanisms for spin relaxation, 
such as the admixing of spin and 2p orbital states, were 
unable to explain the key features of the experimental 
lifetime data, and so introduced a new paradigm that 
includes admixing of valley and spin eigenstates. Our 
theory, which showed good agreement with experiment, 
implies that spin relaxation via phonon emission due to 
spin-orbit coupling can occur in realistic quantum dot 
systems, most likely due to interface disorder. 

Our results show that by tuning of the valley splitting 
in silicon quantum dots, it is possible to always ensure the 
long lifetimes (Ti > 1 s) required for robust spin qubit 
operation. Despite this, the excited valley state will gen- 
erally be lower than orbital states in small quantum dots, 
placing an ultimate limit on the lifetimes accessible in 
very small dots, due to the spin- valley mixing described 
above. 

Electrical manipulation of the valley states is also 
a fundamental requirement to perform coherent valley 
operations. Hence, this work may pave the road for 
the i mplem entation of recently proposed valley-based 
qubit 4^^^. The experimental relaxation rate at the ob- 
served hot-spot was found to be fast {T^^ > 1 kHz) for 
our devices, implying a fast inter- valley relaxation rate, 
so strategies to reduce this would need to be developed. 

Finally, in the context of realizing scalable quantum 
computers, these results allow us to address questions of 
device uniformity and reproducibility with greater opti- 
mism. Indeed, our work suggests that issues related to 
the wide variability of the valley splitting observed in 
silicon nanostructures to date can shift from the elusive 
atomic level (surface roughness, strain, interface disor- 
der) to the more accessible device level, where gate ge- 
ometry and electrostatic confinement can be engineered 
to ensure robust qubit systems. 



VI. METHODS 

The samples fabricated for these experiments are sil- 
icon MOS planar structures. The high purity, near in- 
trinsic, natural isotope silicon substrate has n+ ohmic re- 
gions for source/drain contacts defined via phosphorous 
diffusion. High quality Si02 gate oxide is 8 nm thick and 
is grown by dry oxidation at 800° C. The gates are de- 
fined by electron-beam lithography, Al thermal evapora- 
tion and oxidation. Three layers of AI/AI2O3 are stacked 
and used to selectively form a 2DEG at the Si/Si02 in- 
terface and provide quantum confinement in all the three 
dimensions. 

Measurements are carried out in a dilution refrigera- 
tor with a base temperature T5 ^ 40mK. Flexible coax- 
ial lines fitted with low-temperature low-pass filters con- 
nect the device with the room-temperature electronics. 
In order to reduce pick-up noise, the gates are biased via 
battery powered and opto-isolated voltage sources. The 
SET current is amplified by a room-temperature tran- 
simpedance amplifier and measured via a fast digitizing 
oscilloscope and a lock-in amplifier for the single-shot and 
energy spectrum experiments, respectively. Gate volt- 
age pulses are produced by an arbitrary wave-function 
generator and combined with a DC offset via a room- 
temperature bias-tee. 
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